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Abstract. The single-pion production reactions pp — ► div + , pp — > ripir + and pp — * ppn° were measured at a 
beam momentum of 0.95 GeV/c (T p m 400 MeV) using the short version of the COSY-TOF spectrometer. 
The implementation of a central calorimeter provided particle identification, energy determination and 
neutron detection in addition to time-of-flight and angle measurements. Thus all pion production channels 
were recorded with 1-4 overconstraints. The total and differential cross sections obtained are compared to 
previous data and theoretical calculations. Main emphasis is put on the discussion of the ppiv channel, 
where we obtain angular distributions different from previous experimental results, however, partly in good 
agreement with recent phenomenological and theoretical predictions. In particular we observe very large 
anisotropies for the ty° angular distributions in the kinematical region of small relative proton momenta 
revealing there a dominance of proton spinflip transitions associated with 7r° s- and d-partial waves and 
emphasizing the important role of n° d-waves. 
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1 Introduction 

Single-pion production in the collision between two nucle- 
ons is thought to be the simplest inelastic process between 
two baryons. Therefore it came as a surprise, when first 
near-threshold data on the ppn° channel revealed its 
cross section to be larger than predicted |2|l^|l4^ by nearly 
one order of magnitude. Meanwhile the near-threshold 
data base has been much improved by exclusive measure- 
ments at TRIUMF (TINA and MINA) 0, SATURNE 
(SPES0) and in particular at the cooler storage rings 
CELSIUS (PROMICE/ WAS A) pTISllQ] . COSY (GEM, TOF) 
PTWTTj and IUCF [E], at the latter also with polarized 
beam and target. However, this data base is still far from 
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being complete and often simple observables like angular 
and invariant-mass distributions of the unpolarized cross 
section are still missing. We report in this paper on mea- 
surements of such distributions at a beam momentum of 
0.95 GeV/c (corresponding to T p =397 MeV) for the re- 
actions pp — ► dir + , pp — > ppir and pp — > npir + . For the 
latter we give here the total cross section only, the differ- 
ential distributions will be discussed in a separate paper. 



2 Experiment 

The measurements have been carried out at the Jiilich 
Cooler Synchrotron COSY using the time-of-flight spec- 
trometer TOF at one of its external beam lines. The setup 
of the TOF detector system is displayed in Fig. 1. At the 
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Fig. 1. Sketch of the short COSY-TOF setup used in this 
experiment showing (a) the full detector arrangement with in- 
serts for (b) the start detector region with start wedges and 
hodoscope, (c) the central hodoscope ("Quirl") as the stop de- 
tector for TOF measurements and (d) the central calorimeter. 



entrance of the detector system the beam - collimated to a 
diameter smaller than 2 mm - hits the LH 2 target, which 
has a length of 4 mm, a diameter of 6 mm and 0.9 /j,m 
thick hostaphan foils as entrance and exit windows |13| . 
At a distance of 22 mm downstream of the target the two 
layers of the start detector (each consisting of . 5 mm thick 
scintillators cut into 12 wedge-shaped sectors) were placed 
followed by a two-plane fibre hodoscope (96 x 96 fibers, 
2 mm thick each ) at a distance of 165 mm from target, 
see Fig. lb. Whereas the start detector mainly supplies the 
start times for the time-of-flight (TOF) measurements, the 
fibre hodoscope primarily provides a good angular resolu- 
tion for the detected particle tracks. In its central part the 
TOF-stop detector system consists of the so-called Quirl, 
a 3-layer scintillator system 1081 mm downstream of the 
target shown in Fig. lc and described in detail in Ref. 
|14j - and in its peripheral part of the so-called Ring, also 
a 3-layer scintillator system built in a design analogous 
to the Quirl, however, with inner and outer radii of 560 
and 1540 mm, respectively. Finally behind the Quirl a 
calorimeter (Fig. la,d) was installed for identification of 
charged particles and of neutrons as well as for measuring 
the energy of charged particles. The calorimeter, details 
of which are given in Ref. consists of 84 hexagon- 

shaped scintillator blocks of length 450 mm, which suffices 
to stop deuterons, protons and pions of energies up to 400, 
300 and 160 MeV, respectively. The energy calibration of 
the calorimeter was performed by the detection of cosmic 
muons. 

In the experiment the trigger was set to two hits in 
Quirl and/or Ring associated with two hits in the start 



detector. From straight-line fits to the hit detector ele- 
ments tracks of charged particles are reconstructed. They 
are accepted as good tracks, if they originate in the target 
and have a hit in each detector element the track passes 
through. In this way the angular resolution is better than 
1° both in azimuthal and in polar angles. If there is an 
isolated hit in the calorimeter with no associated hits in 
the preceding detector elements, then this hit qualifies as 
a neutron candidate (further criteria will be discussed be- 
low). In this case the angular resolution of the neutron 
track is given by the size of the hit calorimeter block, i.e. 
by 7 - 8°. By construction of the calorimeter a particle 
will hit one or more calorimeter blocks. The number of 
blocks hit by a particular particle is given by the track 
reconstruction. The total energy deposited by this parti- 
cle in the calorimeter is then just the (calibrated) sum of 
energies deposited in all blocks belonging to the particular 
track. 

In order to have maximum angular coverage by the 
detector elements and to minimize the fraction of charged 
pions decaying in flight before reaching the stop detectors, 
the short version of the TOF spectrometer was used. In 
this way a total polar angle coverage of 3° < €>Lab < 
49° was achieved with the central calorimeter covering the 
region 3° < ® Lab < 28°. For fast particles the 4% energy 
resolution of the calorimeter is superior to that from TOF 
measurements, the resolution of which is reduced by the 
short path length. However, the TOF resolution is still 
much better than the AE resolution of the quirl elements. 
Hence, for particle identification, instead of plotting AE 
versus E ca i, the uncorrected particle energy deposited in 
the calorimeter, we utilize the relation AE ~ (zj (3) 2 with 
the particle charge z = 1 and plot 1 /(3 2 versus E ca i , where 
the particle velocity (3 = v/c is derived from the TOF 
measurement. 

Fig. 2 shows the 1 / 1 — E ca i scatterplot for two-track 
events. The bands for d, p and 7r are well separated. The 
horizontal shadow region on the left of the deuteron band 
stems from deuteron breakup in the calorimeter. Note that 
E ca i in Fig. 2 is not yet corrected for energy and particle 
dependent quenching effects. By applying the quenching 
correction as well as the correction for energy loss in the 
preceding detector elements the kinetic energies of the de- 
tected particles are deduced. 

By identifying and reconstructing the two charged tracks 
of an event the exit channels dir + , npir + and ppir can be 
separated. Kinematically the maximum possible labora- 
tory (lab) polar angles are « 9° for deuterons and » 32° 
for protons (and neutrons). Hence 88% of the angular cov- 
erage for deuterons and 86% of that for protons from single 
pion production are within the angular acceptance of the 
calorimeter. For charged pions the angular coverage has 
been much lower with this setup, since kinematically they 
can extend up to Ohab — 180°. Hence within the angular 
coverage of Quirl and Ring the angular acceptance for 7r + 
has been ~40% only. Nevertheless most of the phase space 
part necessary for a full coverage of the physics in single 
pion production has been covered (see below) by these 
measurements due to the circumstance that the center- 
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Fig. 2. Z\E-E plot for particles stopped in the calorimeter. 
Shown is the first (smaller angle) track of two-track events. For 
the plot the AE information is taken from the TOF information 
obtained by Quirl and start detectors and plotted as 1//3 2 , 
where (3 = v/c denotes the particle velocity normalized to c. 
For the E information the energy deposited in the calorimeter 
E ca i has been taken. 



of-mass (cm) angular distributions have to be symmetric 
about 90° because of identical collision partners in the 
incident channel. 

For 6>Lab < 25.5° the ppir events have been identified 
by requiring the two charged tracks of the events to be 
protons identified in the Quirl-calorimeter system and by 
the condition that their missing mass MM pp meets the 
constraints 100 MeV/c 2 < MM PP < 180 MeV/c 2 (Fig. 
3). For &Lab > 28° the protons of the ppn channel have 
been detected in the Ring. Since kinematically they have 
a small (3 they are easily distinguished from pions hit- 
ting the Ring and also from elastically scattered protons. 
The same applies for the angular region 25.5° > 0Lab < 
28° , where the protons hit only the edge of the calorime- 
ter and in general no longer stop there. In both cases the 
proton energy is calculated from the corresponding TOF 
measured by the Start-Ring and Start-Quirl detector sys- 
tems, respectively. Thus the full kinematically accessible 
angular range was covered for this reaction channel with 
exception of the beam-hole region (0Lab < 3°). 

The dir + events have been selected by identifying both 
deuteron and pion, if both hit the calorimeter or by iden- 
tifying only the deuteron, if the pion hits the Ring. Cor- 
responding missing mass checks have been applied. In ad- 
dition, the coplanarity condition 170° < A<P < 180° is 
used to further distingiush dir + events from three-body 
background. That way even the deuteron identification in 
the calorimeter may be omitted, thus allowing to check 



how well the deuteron breakup in the calorimeter is un- 
der control in the Monte Carlo (MC) simulations. Within 
uncertainties both ways lead to identical results. 

Finally the npir + channel is selected by identifying pro- 
ton and pion in the calorimeter or only the proton in the 
calorimeter, when the second charged track is in the Ring. 
In addition the missing pir mass MM p7T has to meet the 
condition 900 MeV/c 2 < MM p7r < 980 MeV/c 2 . Also 
to suppress background from the dir + channel - in par- 
ticular when the deuteron is broken up and appears as a 
proton in the calorimeter - the pir + track is required to 
be non-coplanar, i.e. A<P < 170° complementary to the 
coplanarity condition given above. Further on the neu- 
tron 4-momentum is reconstructed from the 4-momenta 
of proton and pion and it is checked, whether a calorime- 
ter block in the corresponding (<9, <P) region recorded a 
hit accompanied without any entries recorded in the pre- 
ceding detector elements of the Quirl. If these conditions 
are met, a neutron track is assumed. That way <9„ and <P n 
are determined by the location of this calorimeter block. 
Thus having only the neutron energy undetermined ex- 
perimentally we end up with 3 kinematic overconstraints 
for this channel, whereas we have 1 overconstraint for the 
ppir channel and 4 overconstraints for the dir + channel. 
Corresponding kinematic fits were applied. 

The luminosity of the experiment was determined from 
the analysis of pp elastic scattering. Due to their open- 
ing angle of S pp « 84° between both tracks, such two- 
track events have both hits in the Ring. They are easily 
identified by using in addition the coplanarity constraint 
170° < A<P < 180°. Fig. 4 shows the measured angular dis- 
tribution in comparison with the prediction from the SAID 
database Hi . The data have been efficiency corrected 
by MC simulations of the detector setup. Adjustment of 
the data in absolute height gives the required luminosity 
(L = 5 * 10 33 cm" 2 and dL/dt = 3 * 10 28 cm^s' 1 time 
averaged) and the absolute normalization of the single- 
pion production cross sections, respectively, which will be 
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Fig. 4. Angular distribution of elastic scattering in the cm 
system. The data of elastic scattering proton pairs recorded in 
the ring have been adjusted in absolute height to the SAID 
[Tri] data base (solid line). 

discussed in the following. 



Fig. 5. Pion angular distribution in the cm system for the 
reaction pp — > d% + . The data of this work ( full circles) are 
compared to the SAID |16j data base (solid line). 

we have varied the MC input assuming either pure phase 
space or some reasonable models for the reaction under 
consideration. Statistical uncertainties are negligible com- 
pared to the systematic uncertainties. 



3 Results 

Due to the identity of the collision partners in the entrance 
channel the angular distributions in the overall center-of- 
mass system have to be symmetric about 90 °, i.e. the 
full information about the reaction channels is contained 
already in the intervall 0° < O cm < 90°. Deviations from 
this symmetry in the data indicate systematic uncertain- 
ties in the measurements. Hence we plot - where appropri- 
ate - the full angular range, in order to show the absence 
of major systematic errors present in our measurement. 



Tab. 1: Total cross sections a tot at T p w 400 MeV 
for the reactions pp — » dir + , pp — ► npir + and pp — > ppir 
evaluated in this work and compared to previous measure- 
ments. 





a to t [mb] 




pp — > dn + 


pp — » npir + 


pp — > ppir 


0.74 (2) 


0.47 (2) 


0.100 (7) 


0.72 Q ) 




0.092 (7) b) 






0.070 (7) c) 



■> Ref. [IE] b ) Ref. c) Ref. 



3.1 pp — > dir + 

Absolute and differential cross sections for this reaction 
channel are very well known from previous experiments. 
Hence we use the analysis of our data for this channel 
primarily as a check of the reliability of our measurement 
and data analysis. In Fig. 5 our results for the ir + angular 
distribution are shown in comparison with the prediction 
from the SAID data base ^j. Since in the measurement 
we cover only angles OJ 71 <90° we show in Fig. 5 only 
the appropriate half of the full angular distribution. Note, 
however, that due to the symmetry of the angular distri- 
bution around cosO^f 1 = most of the physically relevant 
phase space part has been covered in this measurement. 
We find good agreement with the SAID data base both 
in absolute magnitude and in the shape of the angular 
distribution. Since the pions of this channel cover the an- 
gular range of both Quirl and Ring, i.e., stem from Quirl- 
Quirl and Quirl- Ring coincidences, the good agreement 
with SAID assures that there are no significant problems 
with correlating the efficiences of Quirl and Ring. 



3.2 pp — > ppn 1 



The evaluated total cross sections for the three chan- 
nels are given in Tab. 1 together with previous results. 
The uncertainties assigned are based on systematics for 
acceptance and efficiency corrections as obtained by vari- 
ation of MC simulations for the detector response, where 



This channel has received increasing attention since first 
measurements in the threshold region jT2] at IUCF and 
later also at CELSIUS [Z| uncovered the total cross section 
to be nearly an order of magnitude larger than predicted 
theoretically [2EEI- Very recent close-to-threshold mea- 
surements at COSY-TOF revealed the experimental total 
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Fig. 6. Spectra for kinetic energies and polar angles (lab sys- 
tem) of detected proton pairs stemming from the pp — > ppir 
reaction. The shaded areas show the corresponding phase space 
distributions for comparison. Dashed and dotted lines give cal- 
culations with the ansatz eq. (4) and the ansatz of Ref. 0, i.e. 
eq. (3), respectively. 

cross sections to be even larger by roughly 50% fUj than 
previously measured. There it was shown that the pp final 
state interaction (FSI) has a very strong influence on the 
reaction dynamics close to threshold with the consequence 
that a substantial part of the cross section is at small lab 
angles, which were missed in IUCF and CELSIUS mea- 
surements near threshold. At higher energies, 320 MeV 
< T p < 400 MeV, where the influence of the pp FSI de- 
creases more and more, the total cross section data mea- 
sured at TRIUMF 0, SATURNE ©, COSY-GEM [Ej 
and CELSIUS (PROMICE/WASA) are in agreement 
with each other with the exception of a 20% discrepancy 
at T p « 400 MeV between CELSIUS and TRIUMF results. 
Our value, see Tab. 1, is in agreement with the CELSIUS 
result. Recent measurements with polarized beam |6J and 
partly also polarized target (Tj added much to the detailed 
knowledge of this reaction from threshold up to T p « 400 
MeV. 

Despite the wealth of experimental information on this 
reaction there remain a number of problems, which are not 
yet sufficiently settled. E.g., the anisotropy of the pion an- 
gular distribution in the overall cm system, characterized 
by the anisotropy parameter 6, which traditionally 0.6 j is 
defined by 

<t(9 c J) ~ l/Z + b*cos 2 O c J, (1) 

where <9£™ denotes the tt° polar angle in the overall cm 
system, shows a big scatter in the results from different 
measurements, see, e.g., Fig. 11 in Ref.fTT] for T p < 400 
MeV and Fig. 6 in Ref.[6J for higher incident energies. Note 
that in Ref.fTT] the Legendre coefficient a 2 for p- waves is 
plotted according to the ansatz 

a{0 c J) ~ 1 + a 2 * (3 cos 2 O c ^ - l)/2, (2) 

In this paper we will use eq. (2) analogously also for fit- 
ting the experimental proton angular distributions a{Op V ^ 
by replacing <9£™ with 0f, p , the proton angle in the pp 
subsystem. The parameters b and a 2 are related by a 2 = 
2b/ (1 + b) ss 26 for b << 1. We will use the quantity a 2 in 
the following discussion of the angular distributions. 



Fig. 7. Dalitz plots for the invariant mass combinations M 2 n0 
versus M 2 and M 2 versus M 2 as obtained from the data 

PP p7T u p7T u 

for the pp — > pp-K reaction. Note that the plots are efficiency 
but not acceptance corrected, hence the tiny deviations from 
the elliptic circumference at the upper corners due to the ex- 
cluded beam-hole region. 



Near-isotropy is found for T p < 400 MeV with a 2 stag- 
gering between -0.1 and +0.1 - with a tendency for nega- 
tive values. The latter would mean that d-wave contribu- 
tions inducing negative a 2 values are already present close 
to threshold || . A clearer trend towards positive a 2 values 
is observed for T p > 400 MeV. 

As we will show below, our results for a 2 are at vari- 
ance with previous results and hence need a more detailed 
consideration. To this end we start the presentation and 
discussion of our results first with energy and angular dis- 
tribution of the protons in the lab system as displayed 
in Fig. 6. In this and in the following one-dimensional fig- 
ures phase space distributions are shown by shaded areas 
for comparison. We see that the data do not deviate vig- 
orously from the phase space distributions, as we would 
expect, e.g., if A excitation would play a dominant role 
in this reaction channel. This is also visible in the experi- 
mental Dalitz plots of versus M 2 p and M 2 ^ versus 
displayed in Fig. 7. The data cover essentially the 
full available elliptic phase space areas and yield distri- 
butions, which are close to flat with just one pronounced 
excursion in the region of the pp FSI. This is reflected also 
in the projections of the Dalitz plot leading to the spectra 
of the invariant masses M p ^o and M pp (Fig. 8). The latter 
exhibits a small spike at the pp threshold due to pp FSI, 
though in total its influence is of minor importance as ex- 
pected from the small amount of s-wave between the two 
protons being available at this energy In the M p7T o 
spectrum we compare our data with the ones from CEL- 
SIUS as given in Ref.|2j. We find agreement between both 
data sets with the exception of the region around M p7r o « 
1100 MeV/c 2 , where we obtain a somewhat larger yield. 

Now we turn again to the angular distributions, which 
are shown in Fig. 9, on the left for the pions denoted by 
their cm polar angles <9£™ and on the right for protons 
in the pp subsystem (Jackson frame) denoted by Of, p , i.e. 
we use the same coordinate system scheme as defined in 
the IUCF publication [T2|. Both distributions are close to 
flat, exhibit, however, a clearly negative anisotropy pa- 



6 



S. Abd El-Samad et al.: Single-Pion Production in pp Collisions at 0.95 GeV/c (I) 



53 4 
> 

I 

-a 

i 1 



jp 4 

u 
> 



V/; 



1880 1900 1920 
M pp [MeV/c 2 ] 



1080 1100 1120 

lYIJMeV/c 2 ] 



Fig. 8. Differential cross sections in dependence of invariant 
masses M pp and M pn o for the pp — > pp7r° reaction. Data of 
this work are shown by full circles, the ones from Ref. [2] by 
open circles and phase space by the shaded area. For the ex- 
planation of dashed and dotted lines see caption of Fig. 6. The 
dash-dotted lines in the M pp distribution show the partial wave 
contributions Ss + Sd and Pp as obtained from eq. (4). 
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Fig. 9. Angular distributions of pions (overall cm system) and 
protons ( pp subsystem, Jackson frame) for the pp — > ppir° 
reaction. Data of this work are shown by full circles, the fit 
to the data with eq. (2) by solid lines, the results of Ref. [2] 
by open circles and the prediction of Ref. j3J by the dash- 
dotted curve. For the explanation of dashed and dotted lines 
see caption of Fig. 6. 



rameter with a 2 = 
for protons. 



-0.12(1) for pions and a 2 = -0.10(1) 



4 Discussion of Results 

The negative anisotropy parameter observed in this ex- 
periment for the 7T° angular distribution comes as a sur- 
prise, since all previous experiments around T p « 400 
MeV gave - or indicated at least - a positive value for 
the pions, the most serious discrepancy being with the 
PROMICE/WASA results E] of a 2 = +0.127(7) (note 
that the b values given in need to be divided by a 
factor of 3, in order to comply with our definition of b in 
eq. 1), since this measurement provided the best statistics 
and phase space coverage of all previous experiments at 
this energy. Although this work reports negative a 2 values 
for T p < 360 MeV, positive values are found for T p > 
400 MeV. 

From their analysis of polarization data the authors 
of Ref. J2| also find some predictions for the unpolarized 
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Fig. 10. Differential cross section in dependence of q, which 
is half of the relative momentum between the proton pair 
(q 2 — Mpp/A — m p ). Data of this work (full circles) are com- 
pared to results from Ref. 9 j (open circles) and to phase space 
(shaded area). For ease of comparison the data from Ref. [|] 
have been normalized to the maximum of our distribution. For 
the explanation of dashed and dotted lines see caption of Fig. 
6. 



angular distributions, though with very large uncertain- 
ties. For protons they get a 2 = —0.34(81) and for pions 
a 2 = 0.17(11). The first one agrees in sign and value with 
our results, however, their value for the pions has an oppo- 
site sign. Nevertheless, since these numbers were obtained 
only indirectly with very large uncertainties, this is not 
a point of major concern - in particular since 7r d-waves 
, which as we demonstrate in this paper are vital for a 
proper understanding of the reaction, are not taken into 
account in Ref. Jl2}. We also note that a recent measure- 
ment of this reaction at T p = 400 MeV was carried out at 
CELSIUS- WASA, too. The results of its analysis also pro- 
vide a negative a 2 parameter for the pion angular distribu- 
tion |17U18j . If we restrict our data to the same proton an- 
gular range covered in the CELSIUS- WAS A experiment, 
then we find full consistency between both results. 

On the theoretical side extensive meson-exchange cal- 
culations as well as partial wave analyses were carried out 
very recently by the Julich theory group |19U20U21| . Their 
calculations partly are in good agreement with the polar- 
ization data. Interestingly their prediction for proton and 
pion angular distributions (shown in Fig. 9 by the dash- 
dotted lines) also results in negative a 2 agreement with 
our results. More advanced calculations based on chiral 
pertubation theory are in progress. First steps in this di- 
rection have already been taken by this group [22 . 

In order to get a better insight into this problem and 
to compare in more detail with the CELSIUS results we 
next consider the angular distributions in dependence of 
q, defined as half of the relative momentum between the 



two protons and given by q 



M 2 

pp 



/4 



m 



v 



In this 



definition q then denotes the momentum of a proton in 
the pp subsystem. The distribution obtained from our 
data is compared in Fig. 10 with the one given in Ref. (Hj. 
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Fig. 11. Pion angular distributions (overall cm system) for 
4 different (/-ranges as indicated for the pp — > ppir° reaction. 
Data from this work are shown by full circles, the fit to the 
data with eq. 2 by solid lines. For the explanation of dashed 
and dotted lines see caption of Fig. 6. 



Whereas good agreement between both results is found 
for q < 150 MeV/c, substantial deviations appear for 
larger q values. In general large q values are associated 
with large opening angles between the protons, which in 
our experiment are fully covered. However, with the setup 
at PROMICE/WASA, where protons were detected only 
for lab angles smaller than 20°, the high-<? region was not 
covered as completely as at TOF and substantial extrapo- 
lations had to be applied in the acceptance correction for 
this part of the CELSIUS data. 



4.1 g-dependence of angular distributions 

In Figs. 11 and 12 we plot pion and proton angular dis- 
tributions for the four different q regions indicated in Fig. 
10. The corresponding anisotropy parameters a 2 obtained 
by fitting eq. (2) to the data are indicated in each of the 
plots. We see that the pion angular distribution is strongly 
anisotropic with a 2 = — 1 for the lowest q range, where 
there are also results from Ref. P] - and in fact they also 
exhibit a very strong anisotropy with a 2 < 0. With in- 
creasing q the pion angular distributions are getting grad- 
ually flatter and for large q the (9£™ distribution gets even 
curved with a 2 > 0. 

The lowest q range (q < 53 MeV/c) is particularily 
interesting. It corresponds to a kinetic energy of the pro- 
tons in the pp subsystem of T pp = q 2 /m p < 3 MeV, a 
constraint, which has been frequently used to have the pp 



subsystem safely selected in a relative 5-wave state. That 
way especially simple configurations in the exit channel 
are selected allowing a deeper insight into basic reaction 
mechanisms of ir° production. In particular, configurations 
are selected that way, where — » AN (I = 1) —* ppir , i.e., 
which have the AN system in relative p-wave in the in- 
termediate state. In fact, the observed exceptionally large 
pion anisotropy of a 2 = —1 equivalent to a pure sin 2 0^" 1 
distribution means that at T p = 400 MeV the ir° produc- 
tion process associated with S-wave protons in the final 
state happens to be a pure proton spinflip process origi- 
nating from the transitions 3 Po — > : SqS and 3 P 2 — ► Sod, 
where capital letters refer to partial waves in the pp sys- 
tem and small letters to partial waves of 7r° relative to the 
pp system (see e.g. Table I in Ref. for a list of con- 
tributing partial waves). In fact, this very special situation 
has been predicted in the phenomenological model of Ref. 
PP, which - based on the CELSIUS- WAS A measurement 
at T p = 310 MeV - gave already a good description of the 
RCNP data [2Bj at T p = 300, 318.5, 345 and 390 MeV in 
the T pp < 3 MeV range. If we use this model ansatz the 
cross section is given by 



a/PS ~ FSI[Al + (2A B k 2 



Bfk* 



cos 



2(9 cml 



q 2 [C 2 



D 2 k 2 



(3) 



where PS and FSI stand for phase space and FSI fac- 
tors ( for the latter see eq. (2) of Ref. |8J). The pion cm 
momentum k and the momentum q, which is half the rel- 
ative momentum between the proton pair, are given here 
in units of the pion mass. Aq and Bq denote pion s- and 
d-wave amplitudes for the transitions 3 Po — Sos and 
3 P 2 -> 1 S a d - as defined in Ref. [Hj with B = -1.2 A - 
whereas C 2 = and D = 0.98 Aq stand for proton P-wave 
contributions of the transitions 1 Sq — > 3 PqS, 3 Pq — > 3 P\p 
and Pi — > 3 Pjp with j = 0, 1, 2, respectively. 

If we use eq. (3) together with the parameters Aq, Bo, 
C 2 and D 2 as determined in Ref. jH], then we obtain the 
dotted curves in Figs. 6 and 8 - 12. The agreement of these 
calculations with our data is striking for the lowest q bin 
(Figs. 11, 12) if renormalized to the same absolute cross 
section of this bin. Also for the other differential distri- 
butions it is partly surprisingly good - with the exception 
of the proton angular distributions, where the sin 2 pp 
ansatz badly fails. However, the trend in the pion angu- 
lar distribution is correctly reproduced. We return to this 
point at the end of this chapter, section (ii). 

For the pion angular distribution in the lowest-g region 
a very strong anisotropy with a 2 < has recently also 
been observed in COSY-ANKE measurements at T p = 800 
MeV [23j. In a very recent work by Niskanen ,25| these 
anisotropies are explained by a strong energy dependence 
of the forward cross section due to interfering pion partial 
waves. 

For the O pp proton angular distribution, unfortunately, 
no data are shown in Ref. [9j to compare with (there the 
proton angular distribution in the overall cm system is 
shown instead, which is slightly different from the one in 
the pp subsystem). In the low-g region, i.e. in the region 
affected most strongly by the pp FSI, we expect the 1 So 
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Fig. 12. Proton angular distributions (pp-subsystem, Jackson 
frame) for 4 different q-ranges as indicated for the pp — » pj>7T° 
reaction. Data of this work are shown by full circles, the fit 
to the data with eq. 2 by solid lines. For the explanation of 
dashed and dotted lines see caption of Fig. 6. 



partial wave between the two protons to dominate. Indeed, 
we find the PP angular distribution to be compatible with 
isotropy within uncertainties. With increasing q the dis- 
tribution gets more and more anisotropic with a 2 < 0. 

The observed trend in the angular distributions is not 
unexpected. In the low-g region the pp system has kine- 
matically the least internal freedom and hence no chance 
to develop much dynamics involving higher partial waves. 
At the same time the pion has kinematically the largest 
freedom within the ppir system with the possibility to in- 
volve dynamically higher partial waves, which then show 
up in appreciable anisotropics of the pion angular distribu- 
tions. At large q the situation is reversed and the pions are 
kinematically bound to low partial waves, i.e. flat angular 
distributions. The observed anisotropics with negative a 2 
values both for proton and pion angular distributions in 
the lower q range point to the importance of proton spin- 
flip transitions in this process, which are associated with 
A excitation in the intermediate state as well as with ir° s- 
and d- waves in the exit channel. Our results are in support 
of the conclusions in Ref. |B] that 7r° d-waves play an im- 
portant role already close to threshold, obviously favored 
by the possibilty of exciting the A in this manner. 



4.2 partial-wave description of data 

Finally we come back to the partial-wave ansatz in Ref. 
[Hj, which is a simplified version of the more comprehensive 
ansatz in Ref. |2j and which we have seen to work quite 
well for the pion angular distributions, however, failing 
badly for the proton distributions. In order to overcome 
these shortcomings we modify eq. (3) slightly (and more 
plausible, if we look at the partial-wave expansions given 
in Ref. [HUE]) by 



a/PS 



FSI[A 2 + (2A B k 2 + B 2 k 4 )cos 2 0^ m ] 
+ q 2 [C 2 +~k 2 D 2 1 (l+D 2 (cos 2 0PP~l/3))}, 



(4) 



where the latter coefficients for the Pp waves are re- 
lated to the one in eq. (3) and in Ref. jHJ by D 2 = 2/3D 2 
in combination with D 2 = —3/2. Note that in this ansatz 
Ps and Pp contributions are taken into account only very 
rudimentarily in the hope that still their most important 
parts are covered by the ansatz. Also for the Sd contribu- 
tion a cos i 0" n term has been neglected due to its small- 
ness and in order to keep the ansatz as close as possible 
to the one of Ref. jH] ■ 

Since the model ansatz of Ref. || works already very 
well for the pion angular distributions, we do not touch the 
pion s- and d-wave parts, i.e. leave the correlation Bq = 
-1.2^0, take A as a general scale parameter to reproduce 
the integral cross section and adjust C 2 , D 2 , D 2 for best 
reproduction of the proton angular and q distributions. As 
a result we obtain 

So = -1.2 Aq, C 2 = 0, D\ = 0.22 A% and D 2 = -0.34 . 

Actually we get a slightly better description having C 2 = 
-0.006. However, since C 2 as a squared quantity should 
not be negative, we set C 2 = as was done also in Ref. 
[Hj. The resulting values mean that aside from Ss and 
Sd contributions the dominant contribution comes from 
Pp configurations, whereas Ps configurations are of minor 
importance. The decomposition of the total cross section 
into interfering Ss and Sd wave contributions on the one 
hand and Pp contributions on the other hand is shown 
in Fig. 8 for the M pp invariant mass distribution. We see 
that the Ss + Sd wave part adds significantly to the total 
cross section and accounts very well for the observed FSI 
effect in the M pp invariant mass distribution. Note that 
Sd contributions, which turn out here to be crucial for 
the understanding of the neagative cuirvature of the pion 
angular distributions, have not been taken into account in 
the analysis of Ref. • 

Though this ansatz is still very simple compared to the 
full partial- wave ansatz as given in Refs. [9lTl2] . it is obvi- 
ously sufficient to provide a near quantitative description 
of the data both for angular distributions and invariant 
masses (dashed lines in Figs. 6, 8 -12). It is not the aim of 
this work to provide a perfect partial- wave fit to the data. 
We rather put here the main emphasis on revealing the 
dominating partial waves in this reaction, which are re- 
sponsible for the main signatures in the (unpolarized) dif- 
ferential observables. Compared to Ref. [H] the only major 
change is the replacement of the D 2 term with its inap- 
propriate sin 2 OP? dependence by the D\ and D 2 terms, 
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which provide a more appropriate angular dependence for 
the protons. This modification, however, is not surpris- 
ing, since the proton angular dependence has actually not 
been tested in Ref. [Hj- From the successful description of 
our data as well of those treated in Ref. [HI we conclude 
that the ansatz eq. (4) provides an amazingly successful 
description of the ( unpolarized) pp —> ppir data from 
threshold up to T p = 400 MeV. 



5 Summary 

It has been demonstrated that by addition of the central 
calorimeter the COSY-TOF setup is capable of providing 
a reliable particle identification on the basis of the Z\E-E 
technique. In this way the different single pion production 
channels were separated. The results for the dir + channel 
agree well with previous results. For the ppn° channel sig- 
nificant deviations from previous investigations were ob- 
tained for angular distributions as well as for invariant 
mass spectra. It has been demonstrated, that pions and 
protons exhibit g-dependent angular distributions. Hence 
a full coverage of the phase space appears to be manda- 
tory for reliable experimental results on this issue. For the 
lower q region angular distributions with negative <Z2 pa- 
rameter dominate pointing to the importance of proton 
spinflip transitions associated with ir° s- and d-waves. In 
particular, we observe a pure sin 2 distribution for T pp < 
3 MeV, which derives from a special combination of the 
spinflip transitions 3 Pq —> Sqs and 3 i-2 — * Sod. Dif- 
ferent from previous experiments this has been the first 
measurement at T p rj 400 MeV covering practically the 
full reaction phase space. The data thus may serve as a 
reliable basis for a comprehensive phase shift analysis of 
this reaction. 
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